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Abstract 

Developmental processes are governed by a diverse suite of signaling pathways employing 
reversible phosphorylation. Recent advances in large-scale phosphoproteomic methodologies have 
made possible the identification and quantification of hundreds to thousands of phosphorylation 
sites from primary tissues. Towards a global characterization of proteomic changes across brain 
development, we present the results of a large-scale quantitative mass spectrometry study 
comparing embryonic, newborn and adult murine brain. Using anti-phosphotyrosine immuno- 
affinity chromatography and strong cation exchange (SCX) chromatography, coupled to 
immobilized metal affinity chromatography (IMAC), we identified and quantified over 1,750 
phosphorylation sites and over 1,300 proteins between three developmental states. Bioinformatic 
analyses highlight functions associated with the identified proteins and phosphoproteins and their 
enrichment at distinct developmental stages. These results serve as a primary reference resource 
and reveal dynamic developmental profiles of proteins and phosphoproteins from the developing 
murine brain. 
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1. Introduction 

Since its emergence two short decades ago, mass spectrometry-based proteomics has 
quickly matured to become not only a major force in discovery science, but also an 
important analytical tool in the testing of hypotheses formed in all manner of biochemical 
and biological disciphnes [1]. Essential to the rapid maturation of proteomics have been 
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remarkable advances in instrumentation, sample preparation methodology, and informatics. 
All of these have enabled proteomics to make major strides toward accomplishing two of its 
main goals: accurate protein identification and quantification at large or comprehensive 
scales [1]. Furthermore, significant progress has also been made on the identification and 
quantification of protein modifications. Dominant in these efforts has been the 
characterization of protein phosphorylation, owing to its recognized regulatory roles, the 
defined mass of its adduct, and what, in hindsight, can be considered relatively straight- 
forward approaches to enrich for phosphorylated peptide analytes [2]. 

Previously we conducted phosphoproteomic analyses of murine brain at embryonic day 16.5 
(E16.5) using SCX [3]; neonatal (PO) brain using SCX-IMAC [4] and post-natal day 21 
(P21) brain using anti-phosphotyrosine peptide immunoprecipitation [5]. These studies 
identified thousands of phosphorylation sites from primary tissue at distinct stages. Notably, 
when we compared our PO SCX-IMAC dataset with a P21 dataset obtained in a near- 
identical fashion [6], we found that one third of the phosphorylation sites in the PO dataset 
were not found in the P21 dataset even though the P21 dataset was more than twice the size 
[4]. These data are consistent with the intuitive hypothesis that phosphoproteomes are highly 
dynamic across development. Herein we document and quantitatively profile the changing 
proteomes and phosphoproteomes between three developmental stages of murine brain: 
E16.5, PO and P21. These results will serve as a reference dataset facilitating the generation 
of many hypotheses relating to vertebrate brain development, ranging from orchestrated 
developmental control of protein cohorts to site-specific developmental regulatory 
mechanisms. 

2. Experimental 

2.1. Tissue Harvesting and Preparation of Tryptic Peptides 

Mice were procured and treated in accordance with an institutionally-approved lACUC 
protocol. Timed-pregnant or appropriately-aged CD-I mice were ordered from Charles 
River Canada (Saint Constant, QC, Canada). Timed-pregnant (E16.5), PO and P21 mice 
were sacrificed and dissected after a brief isoflurane administration. Several brains from 
each stage were pooled and the pooled brains of a given stage were quickly weighed and 
placed on ice immediately prior to lysis. Three P21 brains were pooled, 8 PO brains were 
pooled and 15 E16.5 brains were pooled. For dimethyl-labeling experiments, pooled brain 
tissue from each stage (-750 mg) was dounce-homogenized in a total of 25 mL ice-cold urea 
lysis buffer (8 M urea, 1 10 mM NaCl, 25 mM TRIS pH 8.0, 25 mM NaF, 10 mM Na4P207, 
1 mM Na3V04, 50 mM (3-glycerolphosphate, 10 pg/mL leupeptin, 5 Hg/mL pepstatin A and 
1 mM PMSF). The homogenate was sonicated on ice for six 30-s intervals, with 30-s rests 
on ice between blasts, using a Kontes 50 W sonication microprobe tip at 50% duty output. 
Insoluble cellular debris was removed by centrifugation at 7,000x g for 30 min at 4 °C. DTT 
(5 mM final concentration) was added to 20 mL of the clarified lysate. The lysates were then 
incubated at 65 °C for 30 min. After cooling to room temperature, iodoacetamide was added 
to 12 mM and lysates were incubated for one hour at room temperature in the dark. Samples 
were diluted in 25 mM TRIS pH 8.0 to a urea concentration of 2 M. 750 [ig of sequencing- 
grade of modified trypsin (Promega, Madison, WI, USA) was added to each 10 ml of lysate 
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for overnight digestion at 37 °C. Peptides were acidified with trifluoroacetic acid (TFA) to 
0.4% and insoluble material was spun down by centrifugation at 4,000x g for 20 min. The 
clarified supernatant from each developmental stage was applied to a separate, prewashed 
(100% acetonitrile (MeCN)) and equilibrated (0.1% TFA) tC18 (1 g) solid phase extraction 
Sep-Pak column (Waters, Milford, MA, USA). Peptides were washed with 0.1% TFA and 
eluted in 40% MeCN, 0.1%TFA. Eluates were frozen at -80 °C and lyophilized. For non- 
labeled peptides used in the anti-phosphotyrosine peptide immunoprecipitation experiment 
shown in Figure IC, E16.5 murine brain was lysed in brain complex lysis buffer (BCLB; 25 
mM Tris pH 7.2, 137 mM NaCl, 10% glycerol, 1% Igepal, 25 mM NaF, 10 mM Na4P207, 1 
mM Na3V04, 1 mM PMSF, 10 pg/mL leupeptin and 5 |Jg/mL pepstatin A). The crude 
homogenate was clarified by centrifugation at 15,000x g and at 4 °C. Eleven hand-poured 
7.5%-20% acrylamide gradient, preparative gels containing 6 mg protein each were cut into 
four distinct molecular regions. Similar regions were combined and subjected to in-gel 
tryptic digestion prior to peptide extraction and peptide immunoprecipitations. A detailed 
description of peptide preparation in this manner was described previously [3]. 

2.2. Dimethyl Labeling of Tryptic Peptides 

For the large-scale analysis, 10 mg (weighed) of dried tryptic peptides from each 
developmental stage were separately dissolved in 5.0 mL of 1 M HEPES, pH 7.5. 200 |iL of 
fresh 4% D2-formaldehyde (Cambridge Isotopes Laboratories, Tewksbury, MA, USA) and 
200 |iL 600 mM NaCNBD3 (Cambridge Isotopes Laboratories) were added to both E16.5 
and P21 peptides. 200 |jL of 4% formaldehyde and 200 |jL of 600 mM NaCNBHs were 
added to PO peptides. Reactions were allowed to proceed for 10 min, followed by a second 
200 |jL addition of each reagent and another 10 min incubation. Reactions were quenched by 
adding TFA to 10% and peptides were desalted over tClS (0.5 g) solid phase extraction Sep- 
Pak columns as described above with the washing being 2.5% MeCN, 0.1% FA and the 
elution with 40% MeCN, 0.1% TFA. Peptide eluates were frozen and lyophilized. For the 
small-scale analysis shown in Figure 1 A, ~2 |jg of each mixture was subjected to LC- 
MS/MS analysis. 

2.3. Anti-Phosphotyrosine Peptide Immunoprecipitations 

^H-dimethyl PO peptides, ^H-dimethyl E16.5 peptides, ^H-dimethyl P21 peptides were 
dissolved in 1.4 mL of peptide immunoprecipitation buffer (PIPB; 50 mM MOPS/NaOH pH 
7.2, 10 mM Na2HP04, 50 mM KCl). 0.7 mL (5 mg) of ^H-dimethyl PO peptides was mixed 
separately with 0.7 mL (5 mg) of ^H-dimethyl P21 peptides or 0.7 mL or ^H-dimethyl E16.5 
peptides. Peptides were rocked at 4 °C for 30 min and insoluble remnants were removed by 
centrifugation in a microcentrifuge (15,000x g) at 4 °C for 15 min. 30 pL of a 50% slurry of 
immobilized a-pYlOO (Cell Signaling Technology, Danvers, MA, USA) was added to each 
mixed peptide solution. The solutions were rocked overnight at 4 °C and immune complexes 
were loaded onto a 200 |iL gel-loading tip, pinched to arrest the resin. The resin was washed 
five times with PIPB and twice with water. The original flow through was collected for 
SCX-IMAC. Phosphotyrosine-harboring peptides were eluted with two applications of 40 
|jL of 0.15%TFA. The eluate was desalted on stage tips as previously described [5] and 
dried prior to LC-MS/MS. The non-quantitative E16.5 peptide IP was conducted as 
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previously described [5] and used the monoclonal 4G10 a-phosphophotyrosine antibody 
(Upstate Biotech/Millipore, Billerica, MA, USA). 

2.4. SCX-IMAC 

The flow through from peptide IPs was subjected to SCX-IMAC as described previously 
[2,3] but the SCX was modified to use salt bumps rather than an HPLC as described [7]. 
Briefly, dried peptides were resuspended in SCX buffer A (7 mM KH2PO4 pH 2.65, 30% 
MeCN) and applied to a polysulphoethyl A solid phase extraction column (PolyLC Inc, 
Columbia, MD, USA) pre-washed with 80% MeCN and then with H2O, and then 
equilibrated with SCX buffer A for 30 min. Peptides were eluted fractionally with 6 mL 
each of SCX solvent B (7 mM KH2PO4 pH2.65, 30% MeCN, 350 mM KCl) adjusted with 
SCX solvent A to 0 mM, 10 mM, 25 mM, 40 mM, 60 mM, 90 mM and 150 mM KCl. 
Eluates were collected, frozen at -80 °C, lyophilized, and then desalted on tC18 (0.5 g) solid 
phase extraction Sep-Pak columns as described above and dried. 

2.5. Mass Spectrometry and Data Analysis 

Dried peptides were suspended in 2.5% MeCN, 2.5% formic acid (FA) and loaded for nano- 
scale microcapillary LC -MS/MS in an LTQ-Orbitrap MS (Thermo Electron, Waltham, MA, 
USA) fitted to a Finnigan Nanospray II electrospray ionization source, a Surveyor HPLC 
pump plus, and a Micro AS autosampler (all from Thermo Electron) essentially as described 
[5]. Briefly, after an isocratic loading for 15 min in solvent A (2.5% MeCN, 0.15% FA) 
peptides were separated on an increasing MeCN gradient (2.5%-35%) with 0.15% FA from 
15 to 60 min on a 100 |jm internal diameter, in-house prepared 13 cm long MagicC18 
reverse phase column (5 |jm, 200 A; Michrom Bioresources, Auburn, CA, USA) with a 
needle tip diameter of -4.5 |jm. Peptide measurements were identified in a top- 10, data- 
dependent fashion, using the SEQUEST algorithm (Thermo Electron V26.12) against the 
mouse IPI database (mouse IPI v3.60) in a target-decoy approach [8], allowing for 
phosphorylation of serine, threonine, tyrosine (h-79. 96633 Da), oxidation of methionine 
(h-15. 99429 Da), carbamiodomethylation of cysteine (h-57. 02146 Da), and heavy mass 
addition to N-termini and lysine (h-6.03766). Peptides were required to incorporate a static 
mass addition of N-termini and lysine residues through dimethyl labeling (h-28.0313 Da). 
Quantification of heavy and light peptide pairs was accomplished by interrogating MS 1 fall 
scans and comparing integral values of the distinct isotopic envelopes for each peptide pair 
using Vista-based software [9,10]. MS runs from each developmental comparison were 
pooled by subset — non-phosphopeptides, serine and threonine phosphopeptides, and 
phosphotyrosine peptides — and initially filtered below a 1% false discovery rate using an 
automated linear discriminant analysis as previously described [6] weighted by Xcorr, ACn, 
MS2 ion intensity, missed tryptic cleavages, precursor PPM, and peptide length. Relative 
confidence in phosphorylation site localization was assessed using the Ascore algorithm 
[11]. Note that in some cases the Ascore program adjusted the site of phosphorylation from 
what was designated by SEQU|jEST. 

2.6. Immunoblots and Antibodies 

Whole brains from E16.5, PO or P21 CD-I mice (similarly pooled as described in Section 
2.1) were dounce-homogenized in BCLB (see above). The crude homogenate was clarified 
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by centrifugation at 15,000x g and at 4 °C. Normalized extracts were run on 7.5% or 10% 
(37.5:1 polyacrylamide/bis-acrylamide) SDS-PAGE gels to achieve separation. 
Immunoblotting was performed using 0.2 yim nitrocellulose membranes, blocked with 5% 
dry milk in Tris-buffered saline (TEST) essentially as described previously [12], and 
incubated in primary antibodies (a-CRMP2, a-GSK3P pY216 from Santa Cruz 
Biotechnology, Santa Cruz, CA, USA; a-CRMP2 pT514, a-DCX, a-DCX pS334 and a-|3- 
Tubulin from Cell Signaling Technology) diluted in 1.5% bovine serum albumin (BSA) in 
TEST overnight at 4 °C. HRP-conjugated secondary antibodies (Chemicon/Millipore, 
Billerica, MA, USA) were diluted in TEST. Enhanced chemo-luminescence was performed 
prior to visualization by exposure to Hyblot CL film (Denville Scientific, Metuchen, NJ, 
USA). 

2.7. Bioinformatics 

DAVID [13,14] was utilized to generate enriched Gene Ontology categories, clusters and 
enriched KEGG [15,16] pathways from developmental comparisons utilizing IPI accessions 
as gene identifiers. Prior to conducting bioinformatics analyses, phosphopeptides or non- 
phosphopeptides from each developmental comparison were divided into separate sub-lists. 
Further, to investigate the contribution of quantitative changes to particular pathways and 
functional ontologies, phosphoproteins sub-lists were created from genes in the top and 
bottom 2.5% of each quantitative phosphopeptide dataset. 

3. Results and Discussion 

3.1. Comparative Murine Brain Proteomics Reveals Developmentally-Dynamic Proteomes 
and Phiosphioproteomes 

Previously we found major differences in the identified phosphoserine and 
phosphothreonine sites when comparing SCX-IMAC-enriched phosphopeptides from PO 
and P21 murine brain [4]. To determine if major differences in tyrosine phosphorylation 
profiles could be similarly observed between developmental states, we conducted a large- 
scale anti-phosphotyrosine immuno-affinity enrichment from E16.5 murine brain and 
compared it to results we previously reported for P21 murine brain [5]. We found that of the 
162 phosphotyrosine sites we identified from E16.5 brain only 64 (40%) were also found in 
our study of phosphotyrosine sites from P21 brain which identified 409 sites [5]. These 
results are graphed in Figure 1 and the phosphotyrosine sites identified from E16.5 murine 
brain are provided in Supplementary Table S 1 . Note we provide in Supplementary 
Document S 1 expanded descriptions of technical terms used throughout Supplementary 
Tables. 

While the identification of specific sites of phosphorylation at distinct developmental stages 
is important and helpful, the results would be stronger if measured in a more quantitative 
way. Therefore, we next aimed to monitor the developmental dynamics of the murine brain 
phosphoproteome, as well as the proteome generally, using a quantitative mass spectrometry 
workflow (Figure 1) involving reductive amination to introduce stable-isotope -based mass 
tags at primary amines in tryptic peptides [17]. Prior to conducting large scale comparisons, 
we first examined the general effectiveness of the approach on a smaller scale. E16.5 tryptic 
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peptides were dimethylated in one case using formaldehyde and sodium cyanoborohydride 
without heavy labels, while in the other case the formaldehyde and sodium 
cyanoborohydride contained deuterium atoms. The end result of the heavy chemical reaction 
generated heavy-labeled peptides with six deuterium atoms per primary amine, thereby 
providing a mass tag distinguishable from light counterparts. We made 1:2, 1:4 and 1:8 
(light:heavy) mixtures of the samples prior to subjecting them to LC-MS/MS, peptide 
identification and quantification. In each case roughly 1400 peptides were identified and 
quantified. Box plots of the Log2-transformed data including the means, the first and third 
quartiles, and the standard deviations are presented in Figure 1 . These data show that the 
first and third quartiles for each mixture did not overlap, suggesting this method would 
provide in our hands relatively strong binary quantification capacity at least -/+ three Log2 
orders. 

We therefore proceeded with large-scale analyses using whole brain from three 
developmental stages, E16.5, PO and P21. We used binary comparisons: E16.5 compared 
with PO and PO compared with P2 1 . PO peptides were light-labeled and used as the common 
reference, with E16.5 and P21 peptides heavy-labeled. As outlined in Figure 1, we first 
subjected the mixed peptide sets to anti-phosphotyrosine immunoprecipitations. The 
unbound peptides were subjected to SCX-IMAC. Both sets of affinity purified peptides, as 
well as the SCX-IMAC flow through, were desalted and subjected to LC-MS/MS and data 
analysis as detailed in the experimental section. 

Between the two binary comparisons, over 150 unique phosphotyrosine sites 
(Supplementary Tables S2 and S3) and over 1500 unique phosphoserine and 
phosphothreonine sites (Supplementary Tables S4 and S5) were identified and quantified. 
The identification overlap between the two binary comparisons was roughly 25% for each 
phosphorylation site dataset (Figure 1). This means that relative quantification at each of the 
three developmental stages is possible for 25% of the phosphorylation sites and relative 
quantification is possible at two developmental stages for 75% of the phosphorylation sites. 
These results are consistent with the distinct differences that we observed between 
developmental states in our non-quantitative comparisons discussed above. For protein 
identification and quantification, we used a minimum cutoff of three peptides. Roughly 40% 
of the over 1300 identified proteins were found in both binary comparisons enabling their 
relative quantification at each of the three stages (Figure 1). The quantified proteins and the 
individual non-phosphopeptides used for quantification are provided in Supplementary 
Tables S6-S8. 

3.2. Immunoblotting of Developmentally-lmportant Phosphoproteins Shows Agreement 
with Quantitative lUIS Data 

Although while still relatively few in number, several specific phosphorylation events have 
been shown to be critical for proper vertebrate brain development as elaborated on 
previously [2-5,12,18,19]. As an orthogonal approach by which to evaluate the quantitative 
mass spectrometry data, we used immunoblotting to sample a few proteins and a few 
specific phosphorylation sites known to be critical in brain development. We examined 
Collapsin Response Mediator Protein (CRMP2) and its phosphorylation at T514 (Figure 2). 
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CRMP2 phosphorylation by Glycogen Synthase Kinase-3 (GSK-3) at T514 follows an 
initial priming phosphorylation by Cyclin-dependent Kinase 5 (CDK5) and these 
phosphorylation events play important roles in neuronal migration in response to guidance 
cues [20]. We also examined by immunoblotting the levels of GSK-3 and its activation loop 
phosphorylation at tyrosine 216, the levels of tubulin, and the critical regulator of brain 
development Doublecortin (DCX) and DCX phosphorylation at S339, presumably by CDK5 
[21]. All of the immunoblotting results are presented in Figure 2 and are accompanied by 
quantitative graphs generated by the mass spectrometry data of these same proteins and 
phosphorylation sites. The agreement between the immunoblots and the quantitative mass 
spectrometry is strong. The DCX results also parallel results we obtained and reported 
previously using absolute quantification (AQUA) mass spectrometry [22] as well as 
immublotting [4]. 

3.3. Proteomes and Phosphoproteomes Display Distinct Profiles at Different Stages of 
Brain Development 

Graphing Log2 fold changes of the phosphopeptides and the non-phosphopeptides between 
each developmental comparison, we were surprised to find distinctly different profiles. We 
found that for the non-phosphopeptides the fold change between E16.5 and PO brain formed 
a curve with tighter Gaussian characteristics as compared to the phosphopeptides, whose 
fold changes were more spread and leaned toward higher abundances at PO (Figure 3). In 
contrast, the phosphopeptides in the comparison between PO and P21 showed tight Gaussian 
characteristics with the non-phosphopeptides being more spread and leaning slightly toward 
higher abundances at P21 (Figure 3). These results clearly show the highly dynamic 
proteomes and phosphoproteomes across brain development and surely represent the sum of 
the multitude of regulatory mechanisms occurring at these stages including differential 
transcription, translation, and cell differentiation. The more uniform phosphopeptides 
between PO and P21 may reflect that tissues types are differentiating less between these 
stages and phosphorylation events are therefore more drivers of homeostasis than 
developmental processes. 

3.4. Phosphotyrosine Sites Showing the Most Extreme Differences between E16.5 and P21 
Brain Offer Avenues for Distinct Hypothesis Testing 

Large-scale phosphoproteomic analyses have come of age and now greatly out-pace our 
ability to examine the functional consequences of individual phosphorylation events, not to 
mention cohorts! However, curation of phosphorylation sites, as is done at a multitude of 
sites such as PhosphositePlus [23] and Prosite [24] offers investigators interested in 
conducting functional studies, opportunities to envisage hypotheses of mechanistic 
operation. Table 1 shows the phosphotyrosine sites showing the largest (11 fold or more) 
difference in the relative abundance between E16.5 and PO as well as between PO and P21. 
Each one of these sites has been identified in dozens, and in some cases hundreds, of large- 
scale studies. However, in PhosphositePlus only two are reported to be functionally 
characterized: Dabl pY232 which we found previously to be essential for Dabl's interaction 
with Crk/CrkL in Reelin signaling [18], and GSK3 pY216 which has been extensive 
characterized as the activating event in GSK3's activation loop. The large, observed changes 
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in these phosphotyrosine sites invites investigation into their potential roles governing 
vertebrate brain development. 

3.5. Bioinformatics Reveals Functional and Developmental Dynamics in Brain Proteomes 
and Phosphoproteomes 

To determine if specific functional pathways based on protein counts were overrepresented 
relative to the mouse proteome in either our brain phosphopeptide dataset or our non- 
phosphopeptide dataset we queried the Database for Annotation, Visualization and 
Integrated Discovery (DAVID) to obtain Kyoto Encyclopedia of Genes and Genomes 
(KEGG) [15,16] pathway data. The analysis identified several pathways that were enriched 
in the E16.5:P0 dataset and the P0:P21 datasets. The enriched pathways, unique identifiers 
for the proteins in these pathways, and relevant KEGG statistics are presented in 
Supplementary Tables S9 and SIO. The enriched pathways showed important differences 
comparing developmental stage and when considering non-phosphopeptide and 
phosphopeptide datasets, with some pathways enriched exclusively in only one dataset. A 
mechanism to visualize such presence/absence binary differences as well as differences 
where an enriched pathway is found in both datasets but its relative enrichment competes in 
magnitude with other enriched pathways is a simple rank order analysis. To accomplish this 
we ranked the enriched pathways based on the number of protein counts and visualized 
differences in rank order by drawing lines between the enriched pathway from the E16.5:P0 
dataset to the same enriched pathway in the P0:P21 dataset. The slopes of the lines, as well 
as their color enhancement, accentuate rank order differences (Figures 4 and 5). 

Enrichments offer insight into active pathways, particularly when high protein counts (or 
high rank) are considered. We recognize that a simple rank order discounts various 
weighting options including weighting by protein counts or by fold enrichment. However, 
these approaches can also be conducted using the data in Supplementary Tables S9 and SIO. 
Still, this simple rank order approach is useful in showing variable forms of enrichment in 
our datasets. First, some pathways are only enriched at the phosphopeptide level. For 
example, Mitogen-activated protein kinase (MAPK) signaling and cancer signaling 
pathways were the most dominantly-enriched pathways in both developmental 
phosphopeptide datasets. However, MAPK signaling and cancer signaling pathways were 
not enriched in the non-phosphopeptide datasets. This suggests these proliferation-based 
pathways exhibit significant phospho-regulation at all three of these stages of brain 
development, even if their protein levels do not show enrichment compared to brain 
phosphoproteomes generally. Second, several pathways show enrichment at both the 
phosphopeptide and non-phosphopeptide levels. This includes one of the next most highly- 
ranked pathways, regulators of the cytoskeleton. This pathway also exhibits strong 
enrichment and relatively-high ranking in the non-phosphopeptide datasets. However, axon 
guidance pathways are most strongly enriched in the E16.5:P0 phosphopeptide dataset and 
are only enriched in the P0:P21 non-phosphopeptide dataset. This is not surprising given the 
requirement for axon guidance mechanics during neural development. Finally, those 
enriched pathways showing strong changes in rank order, or enrichment in only one dataset, 
offer insight into bias in developmental dynamics. One example, in addition to the axon 
guidance pathway already discussed, is the spliceosome. Splicing interestingly is only 
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enriched in the E16.5:P0 phosphopeptide dataset, but conversely only enriched in the 
P0:P21 non-phosphopeptide dataset. This suggests strong developmental phospho-regulation 
of the spliceosome, but an increase in total spliceosome protein representation post- 
development, perhaps for routine splicing. Developmentally-regulated splicing has been 
appreciated for some time [25,26]. 

The KEGG analysis described above focused on enrichments found in entire datasets. 
However, given we had acquired quantitative mass spectrometry information, and given our 
primary focus on developmentally-regulated phosphorylation events, we next wanted to 
determine the nature of the phosphoproteins that we had identified and which had a relative 
abundance importantly higher in one developmental stage over another. To do this we took 
our Log2 fold change-sorted datasets described in Figure 3 and reduced them to 
phosphoproteins that were quantitatively in the top or bottom 2.5% of the data. These 
proteins were used to query DAVID to look for clusters of gene ontology categories not 
necessarily exclusively pathways (Supplementary Table SI 1, which includes the relevant 
protein identifications and DAVID statistics). These data for both phosphopeptide datasets 
are shown in Table 2 which abbreviates the data to summarize the major differences 
observed. 

Consistent with developmentally-dynamic phosphoproteomes these DAVID clustering 
results identify differential enrichment in gene ontology categories when analyzing the 
phosphopeptides that exhibit the largest changes in abundance between brain stages. 
Importantly, both the KEGG pathways and other DAVID clustering results are tied to 
distinct proteins and in some cases distinct phosphorylation sites which will facilitate the 
generation of hypotheses that lead to the characterization of proteins with putative roles in 
vertebrate brain development. 

4. Conclusions 

We have conducted a large-scale, quantitative mass spectrometry study of murine brain at 
three developmental stages: E16.5, PO and P21. The relative abundances for thousands of 
phosphorylation sites and proteins across these developmental stages were determined. 
Bioinformatic analyses identified phosphoproteins and proteins from each developmental 
stage with enrichments tied to distinct functionality. Several of these functionalities are 
consistent with nervous system development and development generally. However, much 
remains unexplored and it is anticipated that these results will serve as an important resource 
for developmental biologists as well as for cellular biochemists, particularly those studying 
the roles of phosphorylation on individual proteins or protein cohorts. 

Supplementary Material 

Refer to Web version on PubMed Central for supplementary material. 
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Figure 1. 

Workflow and overview of murine brain developmental proteomic and phosphoproteomic 
comparisons. (A) Reductive amination of tryptic peptides generates mass tags for 
quantitative proteomics. E16.5 tryptic peptides were subjected to reductive amination with 
either formaldehyde (CH2O) and sodium cyanoborohydride (NaBH3CN), or CD2O and 
NaBD3CN, to create light or heavy mass tags respectively. Light to heavy mixtures were 
made as indicated and subjected to LC-MS/MS, peptide identification and quantification. 
Approximately 1,400 peptides were identified and quantified. Box plots show mean, first 
and third quartiles and standard deviations; (B) Quantitative proteomics workflow used 
whole brains from E16.5, PO or P21 mice that underwent urea lysis, trypsinization and 
differential dimethyl labeling as described above. 5 mg heavy-labeled tryptic peptides from 
E16.5 or P21 were each combined with 5 mg light-labeled PO. Mixtures were subjected first 
to anti-pY peptide IPs. Supematants from the IPs were subjected to SCX-IMAC 
phosphopeptide enrichment. The flowthrough from the SCX-IMAC was also retained. The 
peptides bound to each affinity chromatography resin, as well as a portion of the unbound 
peptides were subjected to LC-MS/MS, SEQUEST-based peptide identification. Vista-based 
quantification and Ascore phosphorylation site evaluation as appropriate; (C) Summaries of 
phosphopeptide and protein identifications from our large-scale proteomic comparisons. The 
numbers of phosphotyrosine (pY) and serine/threonine phosphorylation sites (pST) and the 
number of proteins identified from non-phosphopeptides are provided as well as the percent 
of the total identified from each comparison. See text for details. 
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Figure 2. 

Immunoblots parallel quantitative mass spectrometry results. Immunoblots and quantitative 
mass spectrometry were conducted as described in the experimental section, and represent 
the protein and phosphospecific protein levels of the indicated proteins at the indicated 
developmental state. Graphical data on the left represent the fold changes of proteins and 
indicated phosphopeptides between each developmental state as determined by mass 
spectrometry. The sites of phosphorylation are underlined. Errors bars for the mass 
spectrometry data represent the standard deviations of the mean for all non-phosphorylated 
tryptic peptides identified from the indicated protein. Error bars for the phosphopeptides 
represent the standard deviations of the mean from all individual measurements of the 
indicated phosphopeptide. Arrowheads at the top right represent the same molecular 
weights. Molecular weight standards in daltons are indicated. 
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Figure 3. 

Quantitative mass spectrometry results show dynamic proteomic and phosphoproteomic 
profiles at distinct developmental stages of the murine brain. (A) Plot of the E16.5:P0 
datasets; (B) Plot of the P0:P21 datasets. Phosphopeptides and non-phosphopeptides were 
placed into 0.5 Log2 fold-width bins and the numbers in each bin were plotted as a 
smoothened line. 
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Enriched Pathways from 
Protein Phosphopeptides 
Count Rank Order E16.5:P0 

18 MAPK signaling pathway 

17 Pathways in cancer 

16 Regulation of actin cytoskeleton 

IE Axon guidance 

IS Focal adhesion 

13 Chemokine signaling pathway 

12 Spliceosome 

12 Endocytosis 

12 Fc gamma R-mediated phagocytosis 
12 Leukocyte transendotheliai migration 
12 Insulin signaling pathway 
11 Calcium signaling pathway 
11 VEGF signaling pathway 
11 Tight junction 

11 Neurotrophin signaling pathway 
10 ErbB signaling pathway 
10 Alzheimer's disease 
10 Colorectal cancer 

9 Natural killer cell mediated cytotoxicity 

9 T cell receptor signaling pathway 

8 mTOR signaling pathway 

8 Wnt signaling pathway 

8 B cell receptor signaling pathway 

8 Fc epsilon Rl signaling pathway 

8 Long-term poienuanon 

8 Amyotrophic lateral sclerosis (ALS) 

8 Endometrial cancer 

8 Prostate cancer 

7 Phosphatidylinositol signaling system 
7 Adherens junction 
7 Gap junction 

7 Progesterone -mediated oocyte maturation 
7 Type II diabetes mellitus 
7 Glioma 

6 Inositol phosphate metabolism 
6 Apoptosis 
6 Melanoma 

6 NOD-like receptor signaling pathway 

A Aldosterone-regulated sodium reabsorption 

4 Prion diseases 

4 Thyroid cancer 

Not Enriched in E16.5:P0 Comparison 




Enriched Pathways from 
Phosphopeptides 
ink Rank Order P0:P21 

1. MAPK signaling pathway 

2. Regulation of actin cytoskeleton 

3. Pathways in cancer 

4. Calcium signaling pathway 

5. Chemokine signaling pathway 

6. Fc gamma R-mediated phagocytosis 

7. Tight junction 

8. Neurotrophin signaling pathway 

9. Endocytosis 

10. Focal adhesion 

11. ErbB signaling pathway 

12. Leukocyte transendotheliai migration 

13. Axon guidance 

14. Long-term potentiation 

15. Spliceosome 

16. Adherens junction 

17. Insulin signaling pathway 

18. GnRH signaling pathway 

19. Wntsignaling pathway 

20. Vascular smooth muscle contraction 

21. Gap junction 

22. Fc epsilon Rl signaling pathway 

23. Melanogenesis 

24. Phosphatidylinositol signaling system 

25. VEGF signaling pathway 

26. Natural killer cell mediated cytotoxicity 

27. Progesterone-mediated oocyte maturation 

28. Glioma 

29. Oocyte meiosis 

30. Long-term depression 

31. Colorectal cancer 

32. Endometrial cancer 

33. mTOR signaling pathway 

34. Type II diabetes mellitus 

35. Renal celt carcinoma 

36. Prostate cancer 

37. Arrhythmogenic right ventricular cardiomyopathy {, 

38. Non-small cell lung cancer 

39. NOD-like receptor signaling pathway 

40. Aldosterone-regulated sodium reabsorption 

41. Amyotrophic lateral sclerosis (ALS) 

42. Prion diseases 

43. Tliyroid cancer 

44. Not Enriched in P0:P21 Comparison 



Protein 
Count 



Rank Order Change -f 3 to + 6 
Ei6.6:P0 compared to P0:P21 



No Change in Rank Order 
E16.6:P0 compared to P0:P21 



Rank Order Change • 3 to -6 
E16.6:P0 compared to P0;P21 



Rank Order Change + 7 or more 
E16.6:P0 compared to P0:P21 



Rank Order Change +/- 1 to +/- 2 
E16.6:P0 compared to P0:P21 



Rank Order Change - 7 or more 
E16.6:P0 compared to P0;P21 



B Major Rank OrderChanges in Enriched Pathways Based on Phosphopeptides 



Major Decrease in Rank Order 



Major Increase in Rank Order 



(Higher Rank in E16.5:P0 Relative to P0:P21) (Higher Rank in P0:P21 Relative to E16.5:P0) 



Axon guidance 
Spliceosome 
VEGF signaling pathway 
Colorectal cancer 
mTOR signaling pathway 
Amyotrophic lateral sclerosis (ALS) 
Prostate cancer 



Calcium signaling pathway 
Tight junction 

Neurotrophin signaling pathway 
Long-term potentiation 
Adherens junction 
Gap junction 

Phosphatidylinositol signaling system 
Glioma 



Enriched Only in E16.5:PQ Comparison 

Alzheimer's disease 

T cell receptor signaling pathway 

B cell receptor signaling pathway 

inositol phosphate metabolism 

Apoptosis 

Melanoma 



Enriched Only in PQ:P21 Comparison 

GnRH signaling pathway 

Vascular smooth muscle contraction 

Melanogenesis 

Oocyte meiosis 

Long-term depression 

Renal cell carcinoma 

Arrhythmogenic right ventricular cardiomyopathy (ARVC) 



Figure 4. 

(A) Rank order differences in enriched KEGG pathways from proteins identified in 
E16.5:P0 and P0:P21 phosphopeptide datasets indicate potential functional developmental 
dynamics; (B) Major differences (>6 or <6) in rank are summarized. 
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Protein 
Count 



Enriched Pathways 

from Proteins 
Rank Order E16.5:P0 



Enriched Pathways 
from Proteins 
Rank Rank Order P0:P21 



Protein 
Count 



68 Huntington's disease 

64 Alzheimer's disease 

60 Oxidative phosphiorylation 

58 Parl<inson's disease 

42 Ribosome 

35 Calcium signaling pathway 

33 Regulation of actin cytoskeleton 

29 Oocyte meiosis 

27 Cardiac muscle contraction 

26 Glycolysis / Gluconeogenesis 

26 Long-term potentiation 

26 Gap junction 

26 Wnt signaling pathway 

24 Citrate cycle (TCA cycle) 

24 Neurotrophin signaling pathway 

24 Tight junction 

22 GnRH signaling pathway 

21 Melanogenesis 

20 Spliceosome 

19 Proteasome 

19 Fc gamma R-mediated phagocytosis 
19 Vascular smooth muscle contraction 
19 Axon guidance 
17 Pyruvate metabolism 

17 Valine, leucine and isoleucine degradation 

16 Amyotrophic lateral sclerosis (ALS) 

16 Long-term depression 

15 Butanoate metabolism 

14 Propanoate metabolism 

14 Arginine and proline metabolism 

13 Adherens junction 

12 Renal cell carcinoma 

11 Prion diseases 

11 Glioma 

10 Alanine, aspartate and glutamate metabolism 

10 Tryptophan metabolism 

10 Glutathione metabolism 

9 beta-Alanine metabolism 

9 Pentose phosphate pathway 

9 Aldosterone-regulated sodium reabsorption 

9 Fatty acid metabolism 

8 Cysteine and methionine metabolism 

8 Fructose and mannose metabolism 

8 SNARE interactions in vesicular transport 

8 Lysine degradation 

5 Synthesis and degradation of ketone bodies 
5 Glyoxylate and dicarboxylate metabolism 
4 Fatty acid elongation in mitochondria 
4 Valine, leucine and isoleucrne biosynthesis 
Not Enriched in E16.5;P0 Comparison 




1. Ribosome 

2. Huntington's disease 

3. Alzheimer's disease 

4. Spliceosome 

5. Parkinson's disease 

6. Oxidative phosphorylation 

7. Regulation of actin cytoskeleton 

8. Proteasome 

9. Neurotrophin signaling pathway 

10. Tight junction 

11. Glycolysis / Gluconeogenesis 

12. Oocyte meiosis 

13. Fc gamma R-mediated phagocytosis 

14. Gap junction 

,15. Systemic lupus erythematosus 

16. Leukocyte transendothelial migration 

17. Citrate cycle (TCA cycle) 

18. Pyruvate metabolism 

19. Cardiac muscle contraction 

20. Adherens junction 

21. Long-term potentiation 

22. Valine, leucine and isoleucine degradation 
,23. Pentose phosphate pathway 

24. Aminoacyl-tRNA biosynthesis 

25. Propanoate metabolism 

26. Butanoate metabolism 

27. Glutathione metabolism 

28. Arginine and proline metabolism 

29. Prion diseases 

30. Cysteine and methionine metabolism 

31. beta-Alanine metabolism 

32. Terpenoid backbone biosynthesis 

33. Synthesis and degradation of ketone bodies 



50. Not Enriched in P0:P21 Comparison 



Rank Order Change + 3 to + 6 
E16.6:P0 compared to P0:P21 



No Change in Rank Order 
E16.6:P0 compared to P0:P21 



Rank Order Change - 3 to -6 
E16.6:P0 compared to P0:P21 



Rank Order Change + 7 or more 
E16.6:P0 compared to P0:P21 



Rank Order Change +/- 1 to +/- 2 
E16.6:P0 compared to P0:P21 



Rank Order Change - 7 or more 
E16.6:P0 compared to P0:P21 
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B Major Rank Order Changes in Enriched Pathways Based on Non-Phosphopeptides 



Major Decrease in Rank Order 
(Higher Rank in E16.5:P0 Relative to P0:P21) 
Cardiac muscle contraction 
Long-term potentiation 

Enriched Only in E16.5:P0 Comparison 

Calcium signaling pathway 
Wnt signaling patiiway 
GnRH signaling pathway 
Melanogenesis 

Vascular smooth muscle contraction 
Axon guidance 

Amyotrophic lateral sclerosis (ALS) 
Long-term depression 
Renal cell carcinoma 
Glioma 

Alanine, aspartate and glutamate metabolism 

Tryptophan metabolism 

Aldoste rone-regulated sodium reabsorption 

Fatty acid metabolism 

Fructose and mannose metabolism 

SNARE interactions in vesicular transport 

Lysine degradation 

Giyoxylate and dicarboxylate metabolism 
Fatty acid elongation in mitochondria 
Valine, leucine and isoleucine biosynthesis 



Major Increase in Rank Order 
(Higher Rank in P0:P21 Relative to E16.5:P0) 

Spliceosome 
Proteasome 

Fc gamma R-mediated phagocytosis 

Adherens junction 

Pentose phosphate pathway 

Glutathione metabolism 

Cysteine and methionine metabolism 

beta-Alanine metabolism 

Synthesis and degradation of ketone bodies 

Enriched Only in P0:P21 Comparison 

Systemic lupus erythematos 
Leukocyte transendothelial migration 
Aminoacyl-tRNA biosynthesis 
Terpenoid backbone biosynthesis 



Figure 5. 

(A) Rank order differences in enriched KEGG pathways from proteins identified in 
E16.5:P0 and P0:P21 non-phosphopeptide datasets indicate potential, functional 
developmental dynamics; (B) Major differences (>6 or <6) in rank are summarized. 
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Table 2 

DAVID gene ontology summary of phosphoproteins showing major fold changes in the quantitative mass 
spectrometry-based phosphopeptide comparisons. 



Enriched in E16.5 


Enriched in PO 


Cytoskeleton 
Neuronal Projection 
Synaptic Transmission 
Cell Signaling 
Splicing 


Protein Complexes 
Cytoskeleton 

Non-Membrane Bound Organelles 


Enriched Gene Ontology Categories of Phosphoproteins Showing Major Fold Changes in the P0:P21 Phosphopeptide Comparison 


Enriched in PO 


Enriched in P21 


Neurogenesis and Differentiation 
Developmental Growth 
Morphogenesis 
Cell Signaling 
Various Metabolism 


Synaptic Transmission 
Ion Transport 
Cognition 
Cell Signaling 

Non-Membrane Bound Organelles 
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